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All-Band 2G+3G Radial Disc-Cone Antennas: 
Design, Construction and Measurements 

N.I. Yannopoulou, P.E. Zimourtopoulos, E.T. Sarris * 

Antennas Research Group, Austria - Hellas [l, 2] 

EECE Dept, Democritus University of Thrace, Hellas [2 and 3] 

Abstract 

We define as "All-Band 2G+3G" any band that includes all 
frequencies allocated to both 2G and 3G services. We define 
as "Radial Disc-Cone Antenna RDCA" any discone antenna with a 
structure of radial wires. The RDCA was theoretically ana- 
lyzed and software simulated with the purpose of computation- 
ally design a broadband model of it. As an application, a 
broadband RDCA for operation from 800 to 3,000 MHz, which in- 
clude all 2G and 3G frequencies, was designed and an experi- 
mental model was constructed and tested. In order to evaluate 
the agreement between theory and practice, mathematically ex- 
pressed measurement error bounds were computed. 

Introduction Research 

In 1945, Kandoian invented The RDCA was theoretically 

the well-known discone an- analyzed as a group of iden- 
tenna, that is a dipole made tical filamentary V-dipoles 
of a disc above a cone [1] . with unequal arms connected 
In 1953, Nail gave experimen- in parallel. The dipoles re- 
tally two naive relations for cline on equiangular vertical 
the discone dimensions [2]. phi-planes around z-axis to 

In 1987, Rappaport de- form a disconical array, 
signed discones using an Fig.lA shows two such copla- 
N-type connector feed [3]. In nar dipoles conformed with 
1993, Cooke studied a discone the apex angle a. Each 
with a structure of radial V-dipole has a total length L 
wires [4]. In 2005, Kim et equal to the sum of arm 
al. presented a double radial lengths r and s plus the gap 
discone antenna for Ultra 9 between its terminals. 
Wide-Band applications [5] . The simulation was based 

In this short paper we on a suite of developed vis- 
present an All-Band 2G+3G RDCA ual tools which are supported 
fed by an N-type/Female/50- by a fully analyzed, cor- 
Ohm connector. rected and redeveloped edi- 
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tion of the original thin- 
wire computer program by J.H. 
Richmond [6] . 

Two arithmetic criteria 
were adopted for the broad- 
band characterization of a 
model : 

(1) 50-Ohm VSWR lower than 2 

(2) Normalized radiation 
intensity U/Umax lower 
than 3 dB on the 
horizontal plane. 

A visual application pro- 
gram was specifically devel- 
oped to design a broadband 
radial discone with bare 
wires of diameter d embedded 
in free space when the wire 
conductivity, the type of 
feeding connector and the 
frequency band are given. 

The program uses the model 
of a radial discone fed by an 
N-type connector shown in 
Fig. 2. Starting with an ap- 
propriate combination of the 
relations given by [2] -[4] 
the program computes by iter- 
ation in terms of wavelength 
A, the geometric characteris- 
tics r, s, g, a, of the 
broadband model, just when 
the criteria are satisfied. 

Fig. IB shows a Ground 
Plane Antenna GPA that was 
designed for reference and 
consists of equal number of 
cone radials s and a vertical 
monopole with height r. 

As a practical application 
of the broadband design, the 
2G+3G band from 800 to 2,500 



MHz was selected to begin 
with and an experimental ra- 
dial discone of copper wire 
fed by N-type connector was 
built, as shown in Fig. 3. 




Fig . 1 : A - RDCA, B - GPA 



z 




Fig. 2: RDCA - Designed Model 
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In order to demonstrate 
the particular behavior of 
the experimental model, the 
2G+3G band was divided as 
follows : 



2G+3G Sub-Bands 
80Q MHz - 2,500 MHz 



Sub-Band 


MHz 


I 


806 - 960 


II 


1,429 - 1,513 


III 


1,710 - 1,900 


IV 


1,910 - 2,025 


V 


2,110 - 2,170 


VI 


2,400 - 2,499 



Our measurement system 
consists of an EM anechoic 
chamber, a network analyzer, 
a number of support instru- 
ments, a set of standard 
loads of factory accuracy and 
a constructed antenna rota- 
tion mechanism with a built 
hardware control unit of its 
step motor. The combined 
characteristics of system 
parts specify a measurement 
band from 600 to 1300 MHz, 
which overlaps with the 2G/3G 
band. Developed control soft- 
ware synchronizes the system 
and collects data using the 
IEEE-488 protocol. 

A developed general mathe- 
matical method expresses the 
measurement error bounds. An- 
other set of developed soft- 
ware applications processes 
the collected data and com- 
putes the error bounds. 



Results 

The consideration of ra- 
dial discone as an array of 
at least eight 8 V-dipoles 
produces a theta-polarized 
vector radiation pattern with 
magnitude a surface almost by 
revolution around z-axis. So 
the radial discone has indeed 
on the horizontal plane xOy 
the basic properties of a 
vertically polarized almost 
omni-directional antenna, 
that is a fact that encour- 
aged the design of a broad- 
band model by using simula- 
tion . 

The application of the 
broadband criteria to 2G/3G 
band resulted to the design 
of a RDCA with the following 
geometrical characteristics: 



All-Band 2G+3G RDCA 
800 MHz - 3,000 MHz 



Geometry 


Units 


d 


1.5 


[mm] 


r 


44 


[mm] 


g 


6 


[mm] 


s 


125 


[mm] 


a 


60 


[°] 



The RDCA operates from 800 
to 3,000 MHz, which exceeds 
that of 2G+3G band. The ac- 
cordingly constructed experi- 
mental radial discone of 
Fig. 3 should be implied with 
a constructional tolerance of 
±0.5 mm and ±0.5°. 
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The broadband model has a 
directivity from about -0.5 
to 2.9 dBd with slope angle 
between -65° and +58°, but 
the directivity gain on hori- 
zontal plane stays very close 
to the desirable value of 0 
dBi, since it changes from -1 
to +1.7 dBi only. Fig. 4 shows 
that the predicted horizontal 
normalized radiation inten- 
sity remains below 3 dB in- 
deed, while it stays above 0 
dB relative to the reference 
antenna in all 2G+3G sub- 
bands indicated by the verti- 
cal gray strips, when both 
are fed by the same 50-Ohm 
Fig. 3: RDCA experimental model source. 

L/A 

co 0.35 0.70 1.05 1.40 1.75^ 

T3 CD 




Frequency, GHz 



Fig. 4: Predicted radiation intensity on horizontal plane. 
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Fig.5A shows the predicted 
normalized radiation patterns 
in dB at the center of each 
sub-band, which confirms the 
horizontal omni-directional 
radiation properties of the 
broadband model. 

At the center frequency of 
950 MHz of the measurement 
band, the predicted and mea- 
sured radiation intensity on 
the three main cuts of the 
radiation pattern are in good 



agreement, as shown in Fig. 
5B . 

This is made clearer by 
the measurement error bounds 
on a vertical plane as shown 
in Fig. 6. 

Fig. 7 shows that the 50- 
Ohm VSWR predicted results 
for the broadband discone are 
below 2 indeed and almost 
covered by the error bounds 
in the measurement band. 




Fig. 6: Measurement error bounds on a vertical plane 
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Fig. 5: A (Up) Predicted normalized radiation intensity patterns 
at the center of each 2G+3G sub-band - B (Down) Normalized 
radiation intensity pattern at the center of measurements band 
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L/A 




0.6 0.8 1 1.3 2 2.5 3 3.2 

Frequency, GHz 



Fig. 7: Standing wave ratio against frequency or ratio of total 

length to wavelength 



Conclusion 

Prediction and experimen- 
tation in the measurement 
band 600 MHz to 1,300 MHz 
proposes a successfully de- 



signed, constructed, and mea- 
sured Radial Disk Cone An- 
tenna RDCA capable to serve 
All-Band 2G+3G applications 
from 800 MHz to 3,000 MHz. 
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EECE Dept, Democritus University of Thrace, Hellas [2] 

Abstract 



An analytical method was developed to estimate errors in 
quantities depended on full one-port vector network analyzer 
(VNA) measurements using differentials and a complex differ- 
ential error region (DER) was defined. To evaluate the 
method, differences instead of differentials were placed over 
a DER which was then analyzed and compared with another com- 
monly used estimated error. Two real differential error in- 
tervals (DEIs) were defined by the greatest lower and least 
upper bounds of DER projections. To demonstrate the method, a 
typical device under test (DUT) was built and tested against 
frequency. Practically, a DER and its DEIs are solely based 
on manufacturer's data for standard loads and their uncer- 
tainties, measured values and their inaccuracies. 



Introduction 

In full one-port measure- 
ments with a VNA of real charac- 
teristic impedance Z G , a DUT 
with impedance Z has a reflec- 
tion coefficient p defined by 

p= (Z-Z 0 )/(Z + Z 0 ) 

and related to its measured 
value m by the bilinear trans- 
formation 

p = (m-D)/[M(m-D) + R] 

in terms of errors D, M and R 
[1] . This transformation can be 
uniquely determined from given 
distinct p n , n = 1, 2, 3 and re- 
spectively known m k , k = n [2] . 



Research 

We considered p n , m k as the 
elements of given ordered 
triples (A, B, C), (a, b, c), 

solved the resulting system 
and appropriately expressed 
its solution by 

F = X cC(B - A) 

D = X abC(A - B)/F 

M = X c(B - A)/F 

R= [l"l (A - B) (a - b) ]/F 2 

where X and |"l produce two 
more terms from the one 
shown, by rotation of the or- 
dered triple elements. These 
errors were then considered 



THURSDAY 31 OCTOBER 2013 



V4—1 7 



FUNKTECHNIKPLUS # JOURNAL 




A/. I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 



as depended on the indepen- 
dent variables p n , m k . There- 
fore, their differentials 
were expressed in the same 
manner by 

dD = [n (a - b) X (B - C)BC dA 
+ X (b - c) 2 (B - A)(C - A) 

BC da] /F 2 

dM = [X (a - b)(c - a)(B - C) 2 dA 

- n (A - B) X (b - c)da]/F 2 
dR = {X [F + 2(a - b)B(A - C)] 

[(B - C) 2 dA n (a - b) 

- (b - c) 2 da n (A - B) ] }/F 3 

After that, the differential 
of p was expressed by 

dp = [-RdD - (m - D) 2 dM - (m - D)dR 
+ Rdm]/[M(m - D) + R] 2 

and was considered depended, 
through dD, dM and dR, on L = 
7 independent variables and 
their independent differen- 
tials: p n , n = 1, 2, 3 and m k , 
k = n or k = 0 with m 0 = m. 

The developed expressions 
were mechanically verified 
using a developed software 
program for symbolic computa- 
tions . 

Manufacturer's data for 
standard loads used in full- 
one port VNA measurements are 
substituted in p n , and for 
their uncertainties in dp n . 
Since Z 0 is real, the domain 
of each p n is the closed unit 
circle [3]. For |p n | = 0 or 1, 
care must be exercised to re- 
strict its differential value 
onto its domain. The VNA mea- 



surements have specified 
bounded ranges for their mo- 
dulus and argument, so that 
the domain of each m k is a 
bounded circular annular with 
its center at the origin O of 
the complex plane. Measure- 
ment data are substituted in 
m k and manufacturer's data 
for measurement inaccuracy in 
dm k . Uncertainty and inaccu- 
racy data outline domains for 

dp n and dm k . If z = | r | e- 1 ^, stands 
for any of the independent 
variables and dz for its dif- 
ferential then the contribu- 
tion of dz to dp is a summa- 
tion term of the form Wdz, 

with W = | U | e^ v , so that 
Wdz = |U|e j(v + <p) d|r| 

+ |u|e j(v + (p+ TT/2) | r | dtp 

where W is in fact a known 
value of the respective par- 
tial derivative and d|r|, dtp 
are the independent real dif- 
ferentials of the complex dz 
in polar form. Each expres- 
sion Wdz outlines a contour 
for a partial DER around O. 
If z ^ 0, the partial DER is a 
parallelogram with perpendic- 
ular sides d | r | and | r | dcp, 
stretched or contracted by |U| 
and rotated by (V + tp) around 
O. If z = p n = 0, the partial 
DER is a circle with radius 
| U | d | r | . Accordingly, a DER 
is the sum of either L paral- 
lelograms or ( L - 1 ) parallel- 
ograms and 1 circle. DER is 
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then a convex set with con- 
tour either a polygonal line 
with 4L vertices at most, or 
a piecewise curve composed of 
4( L - 1) line segments and 4(L 
- 1) circular arcs at most. 
The greatest lower and least 
upper differential error 
bounds are the end-points of 
DEIs for the real and imagi- 
nary parts of dp and result 
from the projections of DER 
for p on the coordinate axes. 
These conclusions can be ge- 
neralized for any other quan- 
tity directly or indirectly 
depended on all, some or just 
one of the above independent 
variables and their differen- 
tials. Thus, the quantity has 
an L-term DER, where 7>L>1. 
For example, the impedance Z 
of a DUT has the 7-term DER: 

dZ = 2Z 0 dp/(l - p) 2 

Results 

All of the following data 
are specified by manufactur- 
ers of the parts for our mea- 
surement system. This system 
operates from 1 to 1300 MHz 
with 100 Hz PLL stability and 
consists of a type-N Z 0 = 50 Q 
network analyzer, a number of 
support instruments and a set 
of standard loads. The stan- 
dards are: a short circuit A, 
a matching load B with re- 
flection coefficient 0.029 
and an open circuit C with 
reflection coefficient 0.99 
and phase accuracy ±2°. In 
the absence of manufacturer's 



data for A we considered its 
uncertainty equal to that of 
C. So, the following values 
were substituted in the de- 
veloped expressions: 

A = -1, 0 < d | A | <0.01, -180° < 
dtp A < -178° or 178° < dtp A < 180°, 

B = 0, | dB | = 0.029, 

C = 1, -0.01 < d | C | <0, -2° < dtpc 
< + 2 ° 

The annular domain for m k 
of VNA is specified from 0 to 
-70 db in modulus and ±180 
degrees in argument. Measure- 
ments m k result with a deci- 
mal floating point precision 
of 4 digits, for both modulus 
and argument. We consider the 
modulus and argument of dm k 
equal to ±1/2 of the unit in 
the last place of the corre- 
sponding mantissa in modulus 
and argument of m k Conse- 
quently, our system produces 
a DER, either for p or Z, as 
a sum of ( L - 1 ) = 6 parallelo- 
grams and 1 circle, with a 
contour of (4L + 4L) = 48 ver- 
tices at most. 

A suite of developed soft- 
ware applications: (i) con- 

trols the system and collects 
the data in terms of fre- 
quency using the IEEE-488 
protocol, (ii) processes the 
collected data and computes 
the vertices of DER and the 
end-points of its DEIs (iii) 
sketches pictures of DER for 
p and Z in terms of the fre- 
quency steps and make a film 
using them as frames. 
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77.8 86.7 95.5 

Input Resistance R, Ohm 



Fig.l: A typical DER for the impedance Z 




0.2 

Fig. 2: DER for the reflection coefficient p and for its 
associated impedance Z against frequency 
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Fig. 3: Greatest lower and least upper differential error 
bounds for resistance R and reactance X against frequency 



A typical resistor with a 
nominal DC impedance of 50 Q 
±20% tolerance was soldered 
on a type-N base connector 
and enclosed in an aluminium 
box to serve as a simple DUT 
for testing its Z from 2 to 
1289 MHz in 13 MHz steps. The 
center frequency f c = 639 MHz 
was chosen to reveal the de- 
tails of the proposed method 
in Fig.l, where the contour 
of a typical DER for Z is 
outlined with small circles 
as its vertices. This contour 
surrounds that of the 4-terms 
DER due to inaccuracy of mea- 
surements (1) and that of 3- 



terms DER for the uncertainty 
of loads (2). A properly cir- 
cumscribed rectangle of DER 
shows graphically how the 
DEIs for R and X result. The 
commonly used error from the 
matching load only is shown 
as a dotted circle. This is 
in fact a 1-term DER which is 
surrounded from the contour 
of the DER by a factor of 
about 125% to 185% in all di- 
rections. Finally, in the 
same figure, 2 7x2 differences 
AZ resulting from the same 
dp n and dm k , dense enough to 
appear as stripes, were 
placed over DER to compare 
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them with differential dZ 
values. Notably, almost all 
of AZ values are belong to 
DER while the computation 
time for these AZ exceeds 
that for DER by more than one 
order of magnitude. To demon- 
strate the method, a set of 
selected DER frames for p and 
Z are shown in Fig. 2, as 
beads on space curved fila- 
ments against frequency. 



Finally, the computed DEIs 
for R and X are shown in 
Fig. 3 against frequency. 

Conclusion 

The proposed method may be 
efficiently used in the same 
way, to successfully estimate 
errors in any quantity de- 
pended on full one-port vec- 
tor network analyzer measure- 
ments . 
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Abstract 



Since S-parameter measurements without uncertainty cannot 
claim any credibility, the uncertainties in full two-port 
Vector Network Analyzer (VNA) measurements were estimated us- 
ing total complex differentials (Total Differential Errors). 
To express precisely a comparison relation between complex 
differential errors, their differential error regions (DERs) 
were used. To demonstrate the method in the most accurate 
case of a direct zero-length thru, practical results are pre- 
sented for commonly used Z-parameters of a simple, two-port, 
DC resistive T-network, which was built and tested against 
frequency with a VNA measurement system extended by two 
lengthy transmission lines. 



Introduction 

It is well known that in 
full two-port VNA measure- 
ments the S-parameters for a 
two-port Device Under Test 
(DUT) are given in terms of 
their 4 measurements m^, i = 1, 

2, j = 1, 2 by 

~~ £[(^11 - D ) / R] [1 + (m 22 
-D' )M'/R' ] - L(m 21 - X)(m 12 
- X' )/(TT' )}/H (1) 

S 21 = {[l+(m 22 -.D')(M' - L ) /R ' ] 
(m 21 - X)/T}/H (2) 

H = [1 + (m^ - D)M/R] [1 + (m 22 

- D 1 )M ' /R ' ] - LL' (m 21 - X)(m 12 

- X ' ) / ( TT ' ) (3) 



S 22 , S 12 have expressions 

that result from (l)-(2) by 
substituting i, j with j, i 
and D, M, R, L, T, X with D' , 
M', R 1 , T', L', X' and vice- 
versa [1] . These 12 quanti- 
ties have been defined as 
system errors [2] . Stumper 
gave non-generalized expres- 
sions for the partial devia- 
tions of S-parameters due to 
calibration standard uncer- 
tainties, in 2003 [3] . Fur- 

thermore, the developed total 
differential errors for full 
one-port VNA measurements [4] 
are also not generalized in 
the two-port case. To the 
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best of the authors' knowl- 
edge, there are no analytical 
expressions for total differ- 
ential errors in full two- 
port VNA measurements. 



Research 

Since S-parameters are 
functions of 16 complex vari- 
ables, their total differen- 
tial errors were initially 
expressed as 



ds ii = {T T ' (i - MS 11 ) [ R ' + M' (m 22 - D' )] (dm 1:L - dD) 

- RR 1 L(1 - L 1 S 1;L ) [ (m 21 - X) (dm 12 - dX' ) + (m 12 - X' ) (dm 21 - dX) ] 

+ M ' T T 1 [ ( m lrL - D ) ( 1 - MS 1:l ) - RS 1:L ] (dm 22 - dD' ) 

- TT'S 1 i(m 11 - D)[R' + M ' (m 22 - D ' ) ] dM 

+ TT'(m 22 -D')[(m 11 -D)(l - MS 1:L ) - RS 1:L ]dM' 

- ( R ' L( 1 - L'S 1:l ) (m 12 - X' ) (m 21 - X) + 

+ TT'S 11 [R' + M ' (m 22 - D ' ) ] )dR 

- (RL(1 - L'S 11 )(m 12 - X' )(m 21 - X) 

-TT'[(m n -D)(l - MS lrL ) - RS lrL ] )dR' 

- RR' (m 12 - X' ) (m 21 - X) [(1 - L'S 1;L )dL - LS^dL' ] 

+ [(m lrL - D)(1-MS 11 ) - RS 1± ][R' +M'(m 22 - D')] 

(T ' dT+TdT ' ) }/P (4) 

dS 21 = {- MT T'S 21 [R' + M' (m 22 - D' )] (dm 1:L - dD) 

+ RR' LL'S 21 (m 21 - X) (dm 12 - dX' ) 

+ R{T ' [R' + (m 22 -D')(M'-L)] + R ' LL ' S 21 (m 12 - X' )}(dm 21 - dX) 

+ T' (R(m 21 - X) (M' - L) - M'TS 21 [R + M(m 1;L - D)] ) (dm 22 - dD' ) 
-TT'S 21 (m 11 -D)[R' + M ' (m 22 - D ' ) ] dM 
+ T' (m 22 - D' ) (R(m 21 - X) - TS 21 [R + M(m 1;L - D) ] )dM ' 

+ { (m 21 - X) (T ' (m 22 - D ' ) (M ' - L) + R ' [T ' + LL 1 S 21 (m 12 - X ' ) ] ) 
-TT'S 21 [R' + M' (m 22 - D' )]}dR 
+ (R(m 21 - X) [T ' + LL'S 21 (m 12 - X' )] 

- T T'S 21 [R + M(m 1;L - D)] )dR' 

+ R(m 21 - X) [R ' L'S 21 (m 12 - X') - T' (m 22 - D' ) ]dL 
+ RR' LS 21 (m 12 - X' ) (m 21 - X)dL' 

- T'S 21 [R + M(m 1;L - D)] [R* + M' (m 22 - D' )]dT 
+ (R(m 21 - X) [R 1 + (m 22 - D ' ) (M ' - L)] 

- TS 21 [R + M(m 11 - D) ] [R ' + M'(m 22 - D' )])dT'}/P (5) 
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P = TT'[R , +M , (m 22 -D , )][R+H(m 11 -D)] 

- RR'LL' (m 12 - X' )(m 21 - X) (6) 



dS 22 and dS 12 resulted from 

(4), (5) with the mentioned 

substitutions. X, X' errors 
stand for crosstalk measure- 
ments. D, M, R (D 1 , M', R 1 ) er- 
rors are uniquely determined 
in terms of 3 standard loads 
A, B, C (A 1 , B 1 , C') and their 
3 measurements a, b, c (a', b', 
c'), by full one-port VNA mea- 
surements, so the number of 
independent complex variables 
increases from 16 to 22. L, T 
(L', T') errors are accurately 
determined after the replace- 
ment of DUT with a direct 
thru (or approximately, if an 
adapter is used instead) in 
terms of new measurements 
tn, ^21 ( t 22 ' ^- 12 ) anc ^ °f previ- 
ously found quantities. Their 
expressions were appropri- 
ately stated as 



L= [I(ab + ct 11 )C(B-A)]/E (7) 

T= (t 21 -X)[n (A - B) (a - b) ]/ 

(E [IcC(B- A)]) (8) 

E = X(ab + Ct 11 )(B-A) (9) 

where X and |"| produce two more 
terms, from the given one, by 
cyclic rotation of the let- 
ters a, b, c (a 1 , b', c') or A, 
B, C (A ' , B ' , C ' ) . In this way, 
each S-parameter has as total 
differential error dS, a sum 
of 22 differential terms: 

16 due to measurement inaccu- 
racies dm^, dX, dX', d t , da, 

db, dc, da', db', dc' and 6 due 
to standard uncertainties 
given by their manufacturer 
dA, dB, dC, dA', dB 1 , dC'. The 
expressions for dD, dM, dR 
(dD 1 , dM', dR') are known [4]. 
The expressions for the rest 
of differential errors were 
developed as 



dL = {X (B - C)(b - t 1:L )(c - t 1;L ) [ (B - C)(b - a)(c - a)dA 

- (b - c)(B - A)(C - A)da] + [n (A - B)(a - b)] dt^j/E 2 (10) 

dT = {X (t 21 -X)(b- c)(B-C)[(t 11 - c)(b- a) 2 B(A 2 +C 2 ) 

+ (b - t lrL )(c - a) 2 C(A 2 + B 2 ) - 2ABC( b - c)(t 1:L (b + c - 2a) 

- be + a 2 )] [(B - C)(b - a)(c - a)dA 

- (b - c)(B - A)(C - A)da]}/(E 2 [X cC(B - A)] 2 ) 

+ [n (A - B)(a- b)]{[(t 21 -X) Xa(B- C)/E]dt 11 

+ dt 21 - dX}/( E [X cC( B - A) ] ) (11) 
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Each complex differential 
error defines a Differential 
Error Region (DER) on the 
complex plane with projec- 
tions to coordinate axes the 
Differential Error Intervals 
(DEIs) [4]. Obviously, any 



quantity differentiably de- 
pendent on the above vari- 
ables has also a DER. For ex- 
ample, after another correc- 
tion to the given S to Z-pa- 
rameters relations [5], the 
Z-DERs are resulted from 



dZu _ 2Z 0 [ ( 1 S 22 ) 2 dS 11 + (1 S 22 )S 21 dS 12 + (1 S 22 )S 12 dS 21 

+ S 12 S 21 dS 22 ] / [(l - S 1:L )(1 - S 22 ) - S 12 S 21 ] ( 12 ) 

dZ 21 - 2Zg[(l - S 22 ) S 21 dS 11 + S 21 dS 12 + (1 - S 11 )(l - S 22 )dS 21 

+ (l-S 11 )S 21 dS 22 ]/[(l-S 11 )(l-S 22 ) -S 12 S 21 ] 2 (13) 



while dZ 22 , dZ 12 result from 

(12), (13) by application of 

the mentioned substitutions. 

Results 

Six calibration standards, 
in pairs of opposite sex, 
were used and their manufac- 
turers' data were substituted 
in the developed expressions: 

A = -1 = A 1 , 0 < d | A | = d | A 1 | < 
0.01, -180° < dcp A = dcp A , < -178° 

or 178° < dcp A = dcp A , < 180°, 

B = 0 = B' , | dB | = 0.029 = | dB ' | , 
C = 1 = C ' , -0 . 01 < d | C | = d | C ' | 
<0 and -2° < dcp c = dcp c , < +2° . 

The inaccuracy of any VNA 
measurement was conserva- 
tively considered as a sym- 
metric interval defined by 
just 1 unit in the last place 
of the corresponding man- 
tissa, both in modulus and 
argument. Consequently, each 



S-DER is a sum of 20 paral- 
lelograms and 2 circles, with 
a contour of 160 vertices at 
most [4] . 

To demonstrate the method, 
a typical T-network of common 
resistors with nominal DC 
values Z 1 = 24.2 Q, Z 2 = 120 Q 

for the horizontal arms and 
Z 12 =l. 1 £2 for the vertical arm, 

were soldered on type-N base 
connectors of opposite sex 
and enclosed in an aluminium 
box, to form a two-port DUT. 

The VNA measurement system 
was extended by two transmis- 
sion lines of 3.66m and 14m, 
respectively, up to the DUT. 
The DUT was tested from 2 to 
1289 MHz in 13 MHz steps. The 
frequency 1003 MHz was se- 
lected to illustrate the pro- 
posed method for S-DERs shown 
in Fig. 1. 
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To study the total differ- 
ential error, dS was expressed 
as dll + dl, where dll is due to 
the uncertainty of 6 stan- 
dards and dl to the inaccu- 
racy of 16 measurements. The 
contribution of these, con- 
servatively considered mea- 
surement inaccuracies to the 
total differential error is 
as much significant as the 
uncertainties of standard 
loads are. For example, com- 
putations for S 12 over the 

whole measurement band show 
that max|dU| and max|dl| con- 
tribute ~35%-80% and ~25%-70% 
to max|dS 12 |, respectively. 

In addition. Fig. 1 shows how 
the projections of each S-DER 
result its real and imaginary 
DEI. To display the variation 
of S-DER against frequency, a 
number of selected S-DER 
frames are shown in Fig. 2 as 
beads on a space-curved fila- 
ment. It is worth mentioning 
that S 1:l -DER (S 22 -DER) was 

greater than it resulted from 
appropriately organized full 
one-port measurements, as it 
was expected. Finally, the 
computed Z-DEIs are shown in 
Fig. 3, along with their LF 
Z-values . 






Fig . 1 : S-DERs at 1003 MHz 
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1 400 800 1300 1 400 800 1300 

Fig.3:Z-DEIs against frequency 



Conclusion uncertainties in any case 

The proposed method may be where the process equations 
efficiently used to estimate (1), ( 2 ) and (4), (5) can find 

application. 
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